Digital Image Processing 


Image Restoration: 
Noise Removal 
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Image restoration 


• Image restoration is the process of 
recovering the original scene from the 
observed scene which is degraded. 

• Different form enhancement- aim of 
enhancement techniques make images 
visually appealing, whereas restoration 
essentially inverts the degradation- more 
objective. 
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Noise and Images 


The sources of noise in digital 
images arise during image 
acquisition (digitization) and 
transmission 

- Imaging sensors can be 
affected by ambient 
conditions 



- Interference can be added 
to an image during 
transmission 
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Noise Model 


We can consider a noisy image to be modelled as follows: 

g(x, y) = fix, y) + tj(x, y) 


where f(x, y) is the original image pixel, tj(x, y) is the noise term and g(x, y) is the 
resulting noisy pixel 

• If we can estimate the model of the noise in an image, this will 
help us to figure out how to restore the image 

• Noise is assumed to be uncorrelated with pixel values and does 
not depend on spatial coordinates 



Images taken from Gonzalez & Woods, Digital Image Processing (2002) 
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There are many different 
models for the image 
noise term tj(x, y): 

-Gaussian 

• Most common model 


- Rayleigh 

- Erlang 

- Exponential 

- Uniform 

- Impulse 


Salt and pepper noise 
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Noise models 


• Noise is characterized by its probability 
density function. Some of which are: 

• Gaussian, Rayleigh, Uniform, Exponential, 
etc. 

• By far Gaussian is the most popular 
model, because: 

- It is present widely in practice 

- Mathematical ease 
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Filtering to Remove Noise 


We can use spatial filters of different kinds 
to remove different kinds of noise 

The arithmetic mean filter is a very simple 
one and is calculated as follows: 

f(x,y) = — YjS(s,t) 

mn (3S)SS V 

This is implemented as the 
simple smoothing filter 

Blurs the image to remove 
noise 
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Other Means 


There are different kinds of mean filters all of 
which exhibit slightly different behaviour: 

- Geometric Mean 

- Harmonic Mean 

- Contraharmonic Mean 
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Other Means (cont...) 


There are other variants on the mean which 
can give different performance 

Geometric Mean: 


/(*, y ) 


(stfeSxy 


mn 


Achieves similar smoothing to the arithmetic 
mean, but tends to lose less image detail 
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Other Means (cont...) 


Harmonic Mean: 


f(x,y) 


mn 


Z — 


Works well for salt noise, but fails for pepper 
noise 


Also does well for other kinds of noise such 
as Gaussian noise 
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Other Means (cont...) 


Contraharmonic Mean: 

Y,§(s,t) Q+l 


f(x,y) 


(s,t)eS xy 

^8(s 9 f) Q 

(stfeSxy 


Q is the order of the filter and adjusting its 
value changes the filter’s behaviour 

Positive values of Q eliminate pepper noise 

Negative values of Q eliminate salt noise 




Images taken from Gonzalez & Woods, Digital Image Processing (2002) 
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Order Statistics Filters 


Spatial filters that are based on ordering the 
pixel values that make up the 
neighbourhood operated on by the filter 

Useful spatial filters include 

- Median filter 

- Max and min filter 

- Midpoint filter 

-Alphatrimmed mean filter 
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Median Filter: 


Median Filter 


/v 

f(x, y) = median{g(s,t)} 

Excellent at noise removal, without the 
smoothing effects that can occur with other 
smoothing filters 

Particularly good when salt and pepper 
noise is present 
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Max Filter: 

/V 

f(x,y)= max {g0,0} 

(s,t)eS xy 

Min Filter: 

/\ 

f(x,y)= min {#0,0} 

(s,t)eS„ 


Max filter is good for pepper noise and min 
is good for salt noise 
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Midpoint Filter 


Midpoint Filter: 


/(*> 3 ') = ^ 


max {g(s,?))+ min {g0,0} 


(s,t)eS 


*y 


(s,t)eS 


*y 


Good for random Gaussian and uniform 
noise 
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Alpha-Trimmed Mean Filter 


Alpha-Trimmed Mean Filter: 

/(*> y ) = —^—7 

mn-d („ )eS , 

We can delete the d/2 lowest and d/2 highest 
grey levels 

So g r (s, t) represents the remaining mn-d 
pixels 
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Adaptive Filters 


• The filters discussed so far are applied to an 
entire image without any regard for how image 
characteristics vary from one point to another 

• The behaviour of adaptive filters changes 
depending on the characteristics of the image 
inside the filter region 

• We will take a look at the adaptive median filter 
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Adaptive filters for handling noise 


Adapts to some statistical measures on a 
local neighborhood, usually mean and 
variance: 

- g(i,j): the noisy image gray-level value. 

- 2 : noise variance in the image. 

- ° n : local mean in neighborhood. 

jjl l 

- 2 : local variance in neighborhood. 

a L 

vi could be estimated from an uniform area 
in the given image. 
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Adaptive local noise reduction filter 


If noise variance is zero Indicates no 
noise return the observed image. 

If local variance is high compared to the 
noise variance presence of edge or a 
sharp feature-^ return a value close to 
observed grey-value. 

If two variances are equals presence of 
noised return the arithmetic mean 
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4 °< Adaptive local noise reduction filter 


Filter equation: 


a 


f(x, y) = g{x,y) -£ [g(x, y) - m L ] 
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Adaptive Median Filtering 


• The median filter performs relatively well on 
impulse noise as long as the spatial density of the 
impulse noise is not large 

• The adaptive median filter can handle much more 
spatially dense impulse noise, and also performs 
some smoothing for non-impulse noise 

• The key insight in the adaptive median filter is that 

the depending on the 

characteristics of the image 
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Adaptive Median Filtering (cont...) 


Remember that filtering looks at each 
original pixel image in turn and generates 
a new filtered pixel 

First examine the following notation: 

- z min = minimum grey level in S 

- z,„ 0 , = maximum grey level in S v 

- z med = median of grey levels in S xy 
-z xy = grey level at coordinates (x, y) 

- S max =maximum allowed size of S 
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Adaptive Median Filtering (cont...) 


Level A: 

■4 1 Z med Z m in 
— Z mec j z max 

IfAl >0 and A2 < 0, Go to level B 
Else increase the window size 

Level B: 

If window size < S max repeat level A 
Else output z 

R / - 7 - 7 

^xy ^’inin 

B2 —7 — 7 

^xy ^max 

If B1 > 0 and B2 < 0, output z xy 

Else output z med 



Images taken from Gonzalez & Woods, Digital Image Processing (2002) 
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Periodic Noise 


•It arises due to 
electrical or electromagnetic 
interference 

• Gives rise to regular noise 
patterns in an image 

• Frequency domain 

techniques in the Fourier 
domain are most effective 
at removing periodic 
noise 
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Band Reject Filters 


• Removing periodic noise form an image 
involves removing a particular range of 
frequencies from that image 

• The reject !ilters can be used for this purpose 


An ideal band reject filter is given as follows: 
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